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Underwater  imaging  poses  significant  challenges  at 
extended  ranges  when  compared  with  similar  situa¬ 
tions  in  the  air.  Even  in  the  clearest  ocean  waters, 
visibility  range  at  best  is  on  the  order  of  tens  of 
meters,  while  the  visibility  in  air  can  be  on  the  order 
of  miles  under  optimal  conditions.  This  significant  re¬ 
duction  in  range  is  caused  by  the  absorbing  and  scat¬ 
tering  nature  of  the  medium,  which  is  the  water  it¬ 
self,  and  constituents  within,  such  as  particles  of 
various  origins,  including  algal  cells,  detritus,  sedi¬ 
ments,  plankton,  and  even  bubbles  near  the  surface. 
Little  attention  has  been  paid  to  the  effects  of  turbu¬ 
lent  flows  on  imaging  outcomes,  despite  evidence  sug¬ 
gesting  that  the  optical  turbulence  can  be  a  limiting 
factor  under  certain  conditions  [1,2],  Overcoming 
such  challenges  to  increase  both  the  reach  and  the 
resolution  is  of  vital  importance  to  military  and  civil¬ 
ian  applications.  It  is  imperative  to  establish  a  good 
understanding  about  the  limiting  factors  under  dif¬ 
ferent  conditions. 

In  underwater  imaging,  particle  scattering  and 
path  radiance  are  commonly  known  dominant  factors 
in  image  contrast  reduction.  Although  it  is  less  intui¬ 
tive  than  the  star  twinkling  caused  by  air  move¬ 
ments,  the  optical  turbulence  has  been  shown  to  limit 
system  performance  in  natural  waters  [1].  Results 
developed  in  previous  atmosphere  research  [3-6]  are 
used  with  modifications  to  reflect  in-water  optical 
conditions.  Key  steps  and  assumptions  are  also  out¬ 
lined  below  for  convenience  and  assessing  the  limita¬ 
tions  of  the  theory.  The  optical  turbulence  in  the 
ocean  is  primarily  caused  by  the  index-of-refraction 
(lOR)  variation  as  functions  of  the  temperature  and 
salinity.  It  has  been  shown  that  lOR  fluctuations  can 
be  expressed  as  linear  combinations  of  individual  el¬ 
ements,  in  terms  of  both  the  power  spectrum  and  the 
structure  function  [7].  Following  the  Kolmogorov 
model  [3],  for  a  fully  developed  turbulent  flow,  under 
the  inertial  subregime,  2itILq<  k<2itIIq  (k*  is  the 
wavenumber  corresponding  to  eddy  scales;  Lq  and  Zq 
denote  outer  and  inner  scale,  respectively  [4]),  the 
power  spectral  density  of  the  lOR  of  the  ocean  waters 
over  the  imaging  range  (r)  can  be  expressed  in  the 
form  of  [7,8] 


(1) 


where  and  it  reflects  the  3D  optical  tur¬ 

bulence  strength.  is  a  constant  and  is  assumed  to 
be  of  order  of  unity  [8].  e  is  the  kinetic  energy  dissi¬ 
pation  rate,  which  typically  ranges  from 
10"^  to  10"^^  m^  s"^  in  natural  waters,  x  relates  to 
the  dissipation  rate  of  temperature  or  salinity  vari¬ 
ances  [8]  and  ranges  from  10"^  to  10"®  s"^  and 

10""^  to  10"^^  psu^  s"^  (psu,  practical  salinity  unit),  re¬ 
spectively  [9,10].  It  is  apparent  that  the  above  equa¬ 
tion  has  the  usual  Kolmogorov  form  found  in  atmo¬ 
spheric  studies  as  <I>^(K',r)  =  0.033C^(r)/c“^^^  [3-5], 
where  again  the  superscript  K  denotes  Kolmogorov 
spectra.  C\  is  the  structure  constant  of  the  lOR  fluc¬ 
tuations,  which  describes  the  optical  turbulence 
strength  at  a  distance  r  from  the  pupil  plane.  is 
equivalent  of  C^,  by  a  constant,  and  ranges  from  10"® 
to  10^  in  the  ocean.  The  above  scalar  relationship  im¬ 
plies  that  the  turbulence  in  water  is  considered  sta¬ 
tistically  isotropic,  homogenous,  and  wide-sense  sta¬ 
tionary  (WSS)  such  that  spatial  autocorrelation 
function  depends  only  on  relative  positions.  It  is  im¬ 
portant  to  recall  that  not  all  turbulent  flows  can  be 
described  by  the  above  spectra  [4,7]. 

It  is  commonly  known  that  spatial  coherence  func¬ 
tions  between  optical  fields  of  any  two  points  can  be 
used  to  describe  the  irradiance  distribution  of  the 
source  image  or  object  [5,11].  Using  the  Wiener- 
Khinchin  theorem  [5],  the  power  spectrum  (1)  is  re¬ 
lated  to  the  spatial  autocorrelation  of  lOR,  which  it¬ 
self  is  directly  linked  to  the  structure  function  of  lOR. 
Combined,  the  optical  transfer  function  can  be  shown 
as  the  equivalent  to  the  spatial  correlation  function 
on  the  pupil  screen.  For  a  time-varying  correlation 
function  under  WSS  conditions,  its  ensemble  average 
can  be  related  to  the  spatial  phase-structure  function 
such  that  the  optical  transfer  function  (OTF)  of  a 
general  incoherent  object  can  be  expressed  following 
the  approach  by  Fried  [4,6].  The  time-averaged,  or 
long  exposure  (LE)  OTF  of  the  optical  turbulence  in 
underwater  environments  takes  the  form 
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OTFle{/)  =  exp 


(2) 


where  \  is  the  mean  wavelength  (530  nm  for  typical 
underwater  transmissions  [12]).  cf,  is  the  distance  be¬ 
tween  the  pupil  plane  and  the  detector;  f  is  the  spa¬ 
tial  frequency  on  the  pupil  plane  in  units  of  inverse 
length.  To  is  the  so-called  seeing  or  Fried  parameter, 
defined  over  the  propagation  distance  r  as 


To  =  0.185 


477^ 


dzCliz) 


‘f 

•^0 


3/5 


(3) 


and  ^  =  27r/\.  If  we  consider  the  optical  turbulence 
throughout  the  imaging  range  to  be  homogenous  and 
isotropic,  it  can  be  described  by  a  constant  indepen¬ 
dent  of  distance  to  the  pupil  plane,  C^(e)  =  C^,  so  that 
fodzC^(2)  =  rCl,  we  have 


Fig.  1.  (Color  online)  A  sinusoid  (single  frequency)  is  used 
to  illustrate  effects  of  the  path  radiance  on  the  modulation 
transfer.  The  left  side  shows  the  case  of  the  initial  wave 
form  (solid)  with  an  amplitude  of  1  (varies  between  0  and  2) 
and  the  wave  form  after  the  transmission  (dashed)  with  an 
amplitude  of  x  from  the  mean,  where  x  also  equals  the 
modulation  (see  text);  the  normalized  path  radiance  (D) 
is  shown  to  elevate  the  previous  wave  form  by  1/(D+1)  for 
all  spatial  frequencies. 


-^max  -^min 
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l+x-(l-x) 
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Adding  the  normalized  radiance  D  received  by  the 
detector,  assuming  nonsaturation  for  simplicity,  the 
modulation  after  the  effects  of  the  path  radiance  now 
becomes 


To  =  0.185 


477^ 
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,  =  0.185 


0.13277^ 
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(4) 


which  is  a  function  of  the  range  (r).  Rq  will  be  re¬ 
ferred  to  as  the  characteristic  seeing  parameter  and 
denotes  the  seeing  at  the  unit  distance  (1  m).  It  is  im¬ 
portant  to  notice  that  the  underwater  seeing  param¬ 
eter  reduces  over  range  at  a  rate  close  to  the  square 
root  of  r,  implying  fast  rolloff  of  high  spatial  frequen¬ 
cies  at  extended  ranges.  Consequently,  the  time- 
averaged  OTFtur  is  also  a  function  of  the  range  r. 


OTFtur(eA,^)  =  exp 


=  exp(-S„0^M, 


(5) 


where  S„  =  3.44(\/i?o)^^  =  1736/iL3\^^  and  the  angular 
spatial  frequency  The  average  wavelength  is 

used  here,  as  the  phase  shift  caused  by  the  lOR 
variation  due  to  the  temperature  or  salinity  is  not 
primarily  wavelength  dependent. 

Another  factor  that  affects  underwater  imaging  is 
the  path  radiance.  This  effect  can  be  quantified  with 
a  simple  setup  that  produces  generalized  solutions.  If 
we  assume  that  one  transferred  frequency  compo¬ 
nent  (the  sinusoid  in  Fig.  1)  has  amplitude  x  above 
the  mean,  which  is  valued  at  1,  it  is  easy  to  see  that 
between  the  original  amplitude  (1)  and  the  current 
value  (x),  the  modulation  function  is 


^max  ^min  1+X  +  /)  (1  X D) 

M*  = - = - 

^max  ^min  l+X  +  /)  +  1  — X  +  /) 

=  MorigMTFpath-  (7) 

The  result  shows  that  the  path  radiance  affects  all 
frequencies  equally  and  that  the  net  effect  is  a  shift 
of  the  entire  modulation  transfer  curve;  thus  it  does 
not  affect  relative  contributions  between  the  turbu¬ 
lence  and  particles.  Naturally,  at  D  =  0  (no  path  radi¬ 
ance),  Eq.  (7)  shows  MTFpath=l  (MTF  is  the  modula¬ 
tion  transfer  function).  One  can  easily  see  MTFp^th 
— >0  with  a  large  D,  describing  the  saturation  by  am¬ 
bient  lights,  which  limits  all  frequencies.  This  is  in 
general  agreement  with  earlier  results  [13]  but  with¬ 
out  any  specific  illumination  limitation.  As  the  ambi¬ 
ent  light  is  incoherent  in  nature,  the  MTF  can  be  con¬ 
sidered  the  same  as  the  OTF. 

Considering  that  random  phase  changes  of  a  wave- 
front  can  be  described  independently  as  a  thin  screen 
that  exists  only  when  the  turbulence  exists,  the  rest¬ 
ing  or  averaged  OTF  will  be  that  of  particles  only. 
This  assures  linearity  of  system  components,  which 
allows  the  cascading  of  OTFs  in  the  frequency  do¬ 
main  [11].  From  this,  we  can  arrive  at  a  simple  un¬ 
derwater  imaging  equation  that  accounts  for  the  par¬ 
ticle  [12],  the  path  radiance,  and  the  turbulence 
scattering, 

OTF(i/r,r)toui  =  OTF(i/',r)pathOTF(i/f,r)parOTF(^i',r)tur 

1  ■ 
2tt0qi}>  ) 
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exp 
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where  Oq  relates  to  the  mean  scattering,  c  and  a  are 
the  beam  attenuation  and  absorption  coefficients,  re¬ 
spectively.  It  is  worth  pointing  out  that  the  OTEp^r 
can  take  many  different  forms,  depending  on  how  the 
scattering  phase  structure  is  incorporated  [14]. 

The  simple  underwater  imaging  model,  Eq.  (8),  ac¬ 
counts  for  scattering  from  particles  and  the  optical 
turbulence  as  well  as  the  path  radiance  in  the  under¬ 
water  environments.  The  primary  aim  of  the  model  is 
to  determine  the  relative  contributions,  which  is  es¬ 
sential  in  assessing  the  limits  of  conventional  passive 
systems  under  different  underwater  conditions.  The 
model  has  been  applied  successfully  to  explain  dis¬ 
crepancies  between  diver  observations  and  model 
outputs,  particularly  involving  high-frequency  com¬ 
ponents  [15].  Selected  optical  conditions  correspond¬ 
ing  to  the  clear  water  (c  =  0.3  m"^),  strong  turbulent 
environments  (e  =  10"^~  10"®,  ;^=  10"^®— 10"^^  i?o 
=0.002^0.004)  are  used  to  illustrate  the  relative  ef¬ 
fects  of  particle  and  turbulence  scattering  on  imaging 
transmissions,  shown  in  Fig.  2.  We  see  that  the  tur¬ 
bulence  scattering  reduces  imaging  details  rapidly, 
especially  with  increased  path  lengths.  It  has  limited 
effects  on  low-frequency  components  when  compared 
with  particle  scattering,  even  in  clear  waters.  Equa¬ 
tion  (8)  and  Fig.  2  also  help  to  explain  why  Mertens 
reported  no  frequency  higher  than  1  cycle/mrad  ob- 


Fig.  2.  (Color  online)  Comparison  of  relative  contributions 
under  different  conditions:  (T),  the  OTF  contribution  from 
the  turbulence;  (P),  the  particle  scattering  contribution; 
(A),  all  combined  contributions.  Figure  legends  under  cor¬ 
responding  labels  indicate  attenuation  coefficients,  imag¬ 
ing  ranges,  and  seeing  parameters  respectively  (in  m“^,  m 
and  m).  The  single-scattering  albedo  of  all  curves  is  as¬ 
sumed  a  constant  (0.8).  The  last  three  curves  in  the  legend 
are  contributions  from  the  particle,  the  turbulence,  and 
combined  (from  top  to  bottom)  under  the  same  optical  con¬ 
ditions.  No  path  radiance  was  included  (0=0). 


served  in  the  field,  which  puzzled  Duntley  [16],  as  at 
such  high  spatial  frequencies,  the  relative  contrast 
decreases  rather  rapidly  towards  zero.  Lastly,  this 
model  helps  to  explain  the  extreme  turbulence  situa¬ 
tion  observed  by  Gilbert  and  Honey  [1].  If  one  con¬ 
verts  the  standard  United  States  Air  Force  (USAF) 
resolution  chart  [1]  line  pairs  to  spatial  frequencies, 
the  first  blurred  group  (-1)  corresponds  to  650  cycles/ 
rad.  Applying  Eq.  (8)  with  c=0.3,  i?o  =  0.0005,  which 
corresponds  to  a  clear  water,  strong  turbulence  con¬ 
dition  (e  =  10"®,  ;^=10"^^),  one  can  see  that  the  total 
contrast  easily  decreases  to  less  than  2%  within  the 
imaging  range,  which  explains  the  complete  disap¬ 
pearance  as  reported  [1]. 

Equation  (8)  reflects  the  optical  properties  of  the 
medium  under  incoherent  cases.  Since  the  coherent 
cutoff  frequency  is  often  less  [5],  this  can  be  used  only 
as  a  crude  estimate  under  coherent  cases.  Further,  it 
includes  neither  effects  of  the  backscattering  nor 
cases  with  saturations  that  only  further  degrade  the 
MTF.  This  approach  is  not  directly  applicable  to  ac¬ 
tive  systems,  such  as  those  gated  and  modulated,  al¬ 
though  modifications  can  be  made  to  reflect  impacts 
on  shorter  integration  time  to  obtain  system  limita¬ 
tions  similar  to  approach  used  in  [6]. 
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